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We investigate numerically magnetic field generation by thermal convection with 
square periodicity cells in a rotating horizontal layer of electrically-conducting fiuid with 
stress-free electrically perfectly conducting boundaries for Rayleigh numbers in the inter- 
val 5100 < R < 5800. Dynamos of three kinds, apparently not encountered before, are 
Q-i' presented: 1) Steady and time-periodic regimes, where the fiow and magnetic field are 

■ symmetric about a vertical axis. In regimes with this symmetry, the global a— effect is 
0> ! insignificant, and the complex structure of the system of amplitude equations controlling 

weakly nonlinear stability of the system to perturbations with large spatial and temporal 
Q I scales suggests that the perturbations are likely to exhibit uncommon complex patterns 

O ' of behaviour, to be studied in the future work. 2) Periodic in time regimes, where mag- 

■ netic field is always concentrated in the interior of the convective layer, in contrast to 
the behaviour first observed by St Pierre (1993) and often perceived as generic for elec- 
trically infinitely conducting boundaries. 3) A dynamo exhibiting chaotic behaviour of 

^ _ heteroclinic nature, where a sample trajectory enjoys excursions between a periodic mag- 

>■ ' netohydrodynamic regime and rolls. The rolls are amagnetic, but generate magnetic field 

^ ■ kinematically. As a result, magnetic energy falls off almost to zero, while the rolls are 

. approached. 
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1 Introduction 



^ I Nonlinear regimes of magnetic field generation by Rayleigh-Benard convection in elect- 

rically conducting fiuid in a plane layer have not yet been studied in sufficient detail 
beyond the first bifurcations of the trivial steady state. The only notable exception 
is the investigation by Podvigina [20, 22], who traced the sequence of bifurcations and 
attractors for some parameter values on increasing Rayleigh number in the absence of 
magnetic field, and examined the ability of the convective fiows to generate magnetic 
field in kinematic and nonlinear regimes. A reasonably complete catalogue of convective 
hydromagnetic (CHM) attractors in various regions of the parameter space, presenting 
their symmetries and other essential properties, remains so far unavailable. From this 
prospective, numerous studies performed at present, where MHD or CHM simulations are 
carried out for extreme parameter values, are of little interest: any emerging structural 
order in the regimes is ultimately destroyed by chaos in the disguise of turbulence. 
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The goal of the present study is to look at a series of CHM regimes not far from the 
onset of convection. We present dynamos of three kinds, to the best of our knowledge not 
detected before: 

• steady and periodic CHM regimes, symmetric about the vertical axis; 

• dynamos with magnetic field concentrating in the middle of the layer with perfectly 

conducting boundaries; 

• a chaotic CHM regime where magnetic field occasionally switches off almost entirely. 
In the remaining part of the introduction we discuss in what respect in our opinion such 
dynamos are remarkable. 

1.1 The symmetry about a vertical axis and parity invariance 

Our main interest lies in convective regimes of magnetic field generation, which are parity- 
invariant or symmetric about a vertical axis, and stable to perturbations of the same 
periodicity. Laminar dynamos possessing various symmetries are worth to be identified, 
for instance, because symmetries play an important role in the analysis of dynamical sys- 
tems and discovery of symmetric regimes can sometimes serve as a basis for a subsequent 
mathematical investigation providing an insight into different mechanisms of magnetic 
field generation. In this subsection we discuss the fundamental reasons, coercing us to 
focus on CHM regimes with the symmetries mentioned in the title. These reasons are 
related to stability properties of CHM regimes with respect to large-scale perturbations. 

Let us remind the definitions of these symmetries. A CHM regime is parity-invariant, 
if V and h are parity-invariant, and 9 is parity-antiinvariant. A vector field f and a scalar 
field / are parity- invariant, if 

f(-x,t) = -f(x,t); /(-x,t)=/(x,t), 

respectively; they are parity-antiinvariant, if in these equations the signs of the right-hand 
sides are reversed. (The centre of symmetry must lie on the mid-plane; without any loss 
of generality we assume that the origin coincides with the centre of symmetry.) A CHM 
regime is symmetric about the vertical axis, if all the three fields v, h and 6 are. A vector 
field f and a scalar field / are symmetric about a vertical axis, if (again, without any loss 
of generality assuming that the origin lies on the axis of symmetry) 

f2{-Xl, -X2, Xs, t) = -f2{xi,X2, X3, t), 
f3{-Xu -X2: Xs, t) = fs{xi, X2, X3, t), 

and 

/(-Xl, -X2, X3) = /(Xi, X2, X3), 

respectively; they are antisymmetric, if in these equations the signs of the right-hand sides 
are reversed. The two symmetries are compatible with the equations (l)-(4) and with the 
boundary conditions (5)- (7) considered in the present study (sec the next section). For 
periodic regimes, the symmetries can also involve a time shift by a half of the period. 

In this and the following paragraph we summarise the findings of Zheligovsky [30, 
31] , who considered weakly nonhnear stabihty of short-scale CHM regimes to large-scale 
perturbations (involving scales that are much larger than the size of the periodicity cell 
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of the perturbed CHM regime). Generically, large-scale stability is controlled by the 
(global) combined kinetic and magnetic a— effect. Equations governing the evolution 
of weakly nonlinear large-scale perturbations (more precisely, of the leading terms in 
expansion in power series in the scale ratio e of averaged in small scales amplitudes of 
neutral small-scale stability modes, comprising the perturbation) are linear; thus, they 
lack any inherent mechanism for saturation of perturbations. Moreover, they involve no 
other terms but the slow time derivative and the operator of the combined a— effect, and 
the spectrum of the latter is symmetric about the imaginary axis. Consequently, the 
equations generically predict a superexponential growth of the perturbations resulting in 
an ultimate destruction of the perturbed regimes on the time scales 0(£~"''). 

However, if any of the two symmetries under consideration are present, the a— effect is 
insignificant. Apparently (see [3, 5]; note that different quantities are used ihid. and in [30, 
31] to measure the strength of the a— effect - they are based on the physical common sense 
in the first case, and ensue from the multiscale asymptotic analysis without a recourse to 
any empirical relation in the second case), magnetic a— effect can also disappear in highly 
turbulent regimes, but simulations do not demonstrate this conclusively. In the absence 
of significant combined a— effect weakly nonlinear large-scale perturbations are controlled 
by a system of nonlinear PDE's for amphtudes of neutral small-scale stabihty modes. The 
system is mixed: equations for the mean magnetic field are evolutionary, the rest ones are 
not; it can involve cubic nonlinearity (see [31]). In contrast to the case of the presence of 
significant a— effect, a priori the equations do not rule out saturation of perturbations at 
finite energy levels. Instability to large-scale perturbations, if present, develops in MHD 
systems with such symmetries on much larger time scales, 0(£~^). Since in the absence of 
significant a— effect weakly nonlinear large-scale perturbations are controlled by a highly 
complex system of nonlinear PDE's, the perturbations are likely to exhibit a variety of 
uncommon complex patterns of behaviour. 

One can look at the same problem at a different angle. Many convective hydrody- 
namic and hydromagnetic regimes are essentially non-space-periodic. The well-known 
examples include the Kiippers-Lortz [10, 4] and the small-angle (see [21, 23] and refer- 
ences therein) instabilities. Consequently, simulations of convective hydrodynamic and 
MHD regimes in periodic boxes are not physically sound. Moreover, convection in a 
horizontally unconstrained domain can "choose" the preferred wave lengths - e.g. on 
the onset of convective motion the width of convective rolls can be uniquely determined 
[6]. One may try to overcome this difficulty in computations by largely expanding the 
box in horizontal directions. The preferred wave lengths of the convective system, which 
are a priori unknown, are approximated the better, the larger is the box; however, the 
fundamental deficiency of this approach (followed, for instance, in [3, 5]) - in that the 
set of wave lengths involved in computations is uniformly discrete - remains invincible. 
Perhaps, employment of asymptotic multiscale solutions in the style of [9, 30, 31] is at 
present the only alternative method, not suffering from this inadequacy, although it has 
its own limitations. 

We are hence in need of examples of CHM regimes possessing any of the above men- 
tioned symmetries to carry out, in the future work, case studies of their large-scale pertur- 
bations. Podvigina [20, 22] reported instances of hydrodynamic steady and time-periodic 
convective regimes, that were symmetric about a vertical axis and stable to short-scale 
perturbations (hence non-chaotic) of the same period, as the perturbed short-scale CHM 
regime; however, they failed to generate magnetic field with the same symmetry. We have 
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found steady and time-periodic CHM regimes, which are symmetric about a vertical axis. 

1.2 Concentration of magnetic field in the middle of a layer 
with perfectly conducting boundaries 

We study convection in a plane rotating layer with stress-free perfectly conducting bound- 
aries held at constant temperatures. Magnetic field generation by thermal convection with 
square periodicity cells in the nonlinear regime was explored for these boundary condi- 
tions in [14, 25, 11, 7, 26, 3, 5, 20, 22]. Simulations [14, 3, 5] were performed for much 
higher Rayleigh numbers (of the order of 10^ — 10^) than the critical ones, resulting in 
the onset of a turbulence-like behaviour of solutions. In the remaining studies cited above 
the Rayleigh number did not exceed a moderate multiple of the critical value and convec- 
tive flows exhibited laminar behaviour, although instances of chaotic behaviour of MHD 
systems - precursors to turbulence - were also obtained. The case of fast rotation was 
considered in [25, 26]; the Rayleigh number was kept below the critical value for the 
onset of non-magnetic convection in [25]. Magnetic field showed a variety of patterns of 
behaviour. 

In the convective dynamos [25, 11] magnetic field concentrated near the boundaries 
in sublayers of the width of a quarter of the layer width. Magnetic energy attained the 
maximum on the boundary in the run labeled N6 in [14] (see a snapshot of magnetic 
field on fig. 6 ibid.; no detailed information is presented about the spatial distribution of 
magnetic field in other runs). St Pierre [25] suggested a physical mechanism for explana- 
tion of concentration of magnetic field near infinitely electrically conducting boundaries: 
if the fluid motion in the vertical direction is not inhibited, magnetic fleld is advected 
to a boundary and remains "locked" there, unable to exit from the fluid layer outside 
the boundary acting as a screen. Thus, these boundary conditions can be beneflcial for 
magnetic field generation. (In dynamo simulations [27] in a plane convective layer and 
[24] in a rotating spherical fiuid shell, the efficiency of dynamo was not significantly af- 
fected by the boundary conditions for magnetic field, which were varied only on the inner 
core boundary in [24].) In the convective dynamos [20, 22] in a non-rotating layer, strong 
magnetic field usually formed ropes near the boundaries, however, Podvigina [22] found 
chaotic regimes, where magnetic field was residing part of the time near the boundaries 
and part of the time inside the layer. Magnetic field computed in [7] had a form of cigar- 
like structures of vertical orientation, mostly located inside the layer, with small blobs 
of strong field near the horizontal boundaries. In the convective dynamo computations 
[26] for a much faster than in our simulations rotation, magnetic field was organised into 
irregular structures with magnetic energy most time residing inside the fiuid layer. 

We have encountered CHM time-periodic regimes, in which, in contrast with all the 
simulations mentioned above, the generated magnetic field always remains concentrated 
in the interior of the convective layer. 

1.3 A heteroclinic CHM regime 

with occasional disappearance of magnetic field 

We have found a chaotic CHM regime of a heteroclinic nature, where a trajectory jumps 
between a CHM periodic orbit and amagnetic steady rolls. Both, the orbit and the 
rolls, are shghtly unstable. The rolls are capable of kinematic magnetic field generation. 
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however, while the trajectory approaches the rolls after the regime has deviated from the 
periodic oscillations, magnetic energy decays almost to zero before the growing magnetic 
mode becomes substantial. 

Prom the viewpoint of the theory of dynamical systems, nothing is particularly sur- 
prising in such a behaviour: heteroclinic cycles between a periodic orbit and two steady 
states were studied in [12] and were observed, for instance, in [13] in a low-order dy- 
namical system, modelling excursions and reversals of the Earth's magnetic field. To the 
best of our knowledge, such behaviour was never before found in numerical simulations 
of "real" physical systems. It is notable, that this CHM regime and regimes observed for 
some sets of parameter values in the VKS (von Karman Sodium) experiment [16] exhibit 
certain similar features, such as intermittency with bursts and extinction of magnetic field 
(see fig. 25 ibid.) and heteroclinic connections of unstable MHD steady states (but the 
regimes of heteroclinic nature apparently do not involve connections to periodic orbits). 



2 Numerical simulation 



Evolution of the Boussinesq CHM system satisfies the non-dimensional Navier-Stokes, 
magnetic induction and heat transfer equations: 



9v 

_ = PV V -Fvxa;-hx(Vxh)-F Prv x 63 + PROea - Vp, 
— = — V^h + Vx vxh, 

at Pm 



de 

dt 



- (v • v)e + v^, 



and the solenoidality conditions 



V-v = 0, 



V 0. 



(1) 
(2) 
(3) 

(4) 



Here h and a; = V x v denote velocity and vorticity, respectively, of a flow of an elec- 
trically conducting fluid, p modified pressure, h magnetic field, 9 = T — [x^ -\- L/2) the 
difference between the normalised temperature T and the linear temperature profile in 
the layer in the absence of fiuid motion; t is time, P, Pm, R and are kinematic Prandtl, 
magnetic Prandtl, Rayleigh and Taylor numbers, respectively, t/2 is the normalised angu- 
lar rotation rate, the unit vector along the axis Xk of the Cartesian coordinate system 
co-rotating with the layer; L is the width of the layer. The CHM system is free: no 
external sources are present in equations (l)-(3). Hence, the regimes are spatially and 
temporally invariant - there is no explicit dependence of any term in the equations either 
on the position in space, or on time. 

We assume that the horizontal boundaries of the layer are stress-free: 



dvi 



dV2 



^3=±L/2 

infinitely electrically conducting: 
dhi 



dx^ 



X3=±L/2 



dx:> 



dh2 

dXr> 



0, 



X3=±L/2 



X3=±L/2 



= 0, 



(5) 



= 0, h. 



X3=±L/2 



^\x3=±L/2 



0, 



(6) 
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and kept at fixed temperatures: 

^ xz=±Ll2 = 0. (7) 

These boundary conditions are the most convenient ones for simulation of Boussinesq 
convection in the layer. They were used in all numerical dynamo simulations cited in the 
Introduction, (5) was assumed in the investigation [21, 23] (also see references therein) of 
the hydrodynamic small-angle instability, and stability regimes to large-scale perturba- 
tions were studied in [30, 31] also for (5)- (7). 

We have carried out a numerical search for examples of symmetric CHM attractors 
with the same period Lh along the horizontal coordinate axes xi and X2- To avoid turbu- 
lisation of the flow and symmetry breaking, we have not used extreme parameter values 
like in some studies cited in the Introduction. Simulations have been performed for L = 1, 
Lh = 10/7, P = 1, r = 91 (i.e. the Taylor number is 8281) for R increasing from 5100 
to 5800. For these values, the critical Rayleigh number for the onset of convective flow in 
the absence of magnetic field is equal to ~5514; for the periodicity box, which is \/2 times 
smaller in size and aligned with the diagonals of the computational periodicity box, the 
critical Rayleigh number is ~5072 (the two critical numbers have been evaluated using a 
precise formula for R derived in [6] for the boundary conditions assumed here). 

As discussed in the Introduction, our prime goal is to find symmetric CHM regimes. 
In the primary (hydrodynamic) bifurcation convective fiows in the form of rolls set in, 
which possess the desirable symmetry. Because of rotation, they can (at least in principle) 
kinematically generate magnetic field, if magnetic Prandtl number is high enough; if the 
dominant magnetic mode also has the symmetry, then the CHM regime emerging in the 
secondary (magnetic) bifurcation remains symmetric. Consequently, we have chosen to 
perform simulations for a (relatively) high value = 100 (in contrast with < 1 used 
for numerical simulations of the Earth's dynamo, characterised by the values ~ 10~^). 

Unhke in [20, 22], we have not attempted to analyse the sequence of occurring bifur- 
cations in detail, or to locate them to a good accuracy, nor to identify all attractors of 
the CHM system for any considered R. 

For the boundary conditions under consideration, the fields are expanded in the Fourier 
series 

Un,i cosn37r(a;3/L + 1/2) 
v= ^ Wn,2 cosn37r(a;3/L + 1/2) e^^ 

n: n3>0 [ Vn,3 Sm risTl {xs/ L + 1/2) 



Sp-(niXi+n2X2) 



(8) 



n: n3>0 



h^ i cosn-^n^x^/ L + 1/2) 
hn,2Cosn3n{x3/L + 1/2) 
hn,3 sin n^irixs/L + 1/2) 



{nixi+n2X2) 



e 



E 



'¥^{niX]_+n2X2) 



9 ^sm n^-K {x3 / L + 1/2) e' 

n: n3>0 

Solenoidality conditions (4) for the series (8) and (9) take the forms 

27rini 27rin2 irn^ 

-Vr,,i H ; — v„2 + ^—'^11,3 = 0; 



(9) 



(10) 



27rini 



27iin2 , 
,1 H ; — hr 



Lh ' Lh ' L 
Standard pseudospectral methods (see [2, 1, 17]) have been apphed with the resolution of 
64 X 64 X 32 Fourier harmonics (after deahasing in computation of advective terms was 
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Table. Regimes of thermal hydromagnetic convection. E^^^, E^^, E'^^: kinetic en- 
ergy per periodicity cell, minimal, average over a period and maximal, respectively; 
-^min> -^max- Same for magnetic energy. SVA: symmetry about the vertical 

axis (+: v and h are both symmetric; — : v is symmetric and h is antisymmetric; 0: v 
is symmetric, no magnetic field generation; No: v and h do not possess the symmetry 
or antisymmetry). P: parity related symmetries (see the text). Type: the type of regime 
(SS: a steady state, DR: steady rolls aligned with the diagonal, with ni=n2 and of the 
same parity in (8), HC: heteroclinic chaos, a number: temporal frequency of a periodic 
orbit) . 
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5100 


1.0 


1.0 


1.0 


0.012 


0.012 


0.012 
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SS 


5140 


2.4 


2.4 
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0.025 


0.025 


0.025 
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SS 


5160 


4.9 
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1.30 


1.46 
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0.18 


5180 
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1.58 
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2.28 
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5190 
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2.10 


2.68 






0.19 
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5.2 


6.8 


16.1 


0.33 


1.86 


12.8 
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6.5 
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6.5 
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8.8 
















DR 


5300 


10.8 


10.8 


10.8 
















DR 


5350 


9.8 


10.0 


10.4 


0.30 


0.55 


0.95 
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0.315 


5400 


11.5 


12.0 


13.0 


0.26 


0.75 


1.85 
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o-i 


0.31 


5500 


15.1 


16.0 


19.0 


0.27 


1.13 


3.71 
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o-i 


0.32 


5550 


16.9 


18.0 


22.9 


0.25 


1.15 


4.68 
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0"! 


0.29 


5600 


20.3 


22.3 


25.5 


1.56 


2.30 


3.77 


No 


0-2 


0.66 


5700 


24.6 


27.3 


31.0 


2.20 


3.20 


4.91 




0"2 


0.76 


5750 


26.3 


29.8 


33.9 


2.28 


3.71 


5.62 


No 


0-2 


0.38 


5800 


26.1 


33.0 


39.7 





3.70 


27.12 


No 


0-3 


HC 



performed). The energy spectra of the computed attractors decay by several (at least 3-4 
for magnetic field, which has the least decreasing energy spectrum) orders of magnitude. 
No significant discrepancies with the results of several test runs with the 128 x 128 x 64 
resolution have been found. 

3 Results of computation 

Properties of the observed CHM attractors are summarised in the table. In the interval 
5210 < i? < 5300 they degenerate into a system of two-dimensional rolls aligned with the 
diagonal of the periodicity box, in which the fiow is independent from the position along 
the axes of the rolls. Although the flows in the rolls has all the three spatial components 
and the Zeldovich antidynamo theorem [28] is not applicable, the rolls do not generate 
magnetic field. In all other saturated regimes that have been simulated, generation takes 
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place. 

For most attractors found in computations, the sums (8)- (10) turn out to be comprised 
of terms, where for some i, j = 1,2,3 the sum rii + rij has a certain parity (see the table). 
Suppose Vn,i 7^ only, if rii + rij is even. This is consistent with the action of the Lorentz 
force (quadratic in h) in the Navier-Stokes equation (1) only, if for all non-zero terms 
in (8) parity of rii + rij is the same - either odd, or even. If a scalar or vector field is 
represented by a sum of the form of (10) or (8), respectively, where rii + n2 is even, then 
the size of its elementary periodicity cell is and the cell sides are oriented along 

the diagonals of the computational periodicity cell. In particular, for R — 5100, 5140 
and 5350 < R < 5700 v^^i ^ only, if ni + n2 is even, and h^^i ^ only, if rii + n2 is 
odd. Thus, the periodicity cell of the fiow is times less than the periodicity cell of 
magnetic field, i.e. a mild scale separation is present in these CHM regimes. If coefficients 
of a field do not vanish only, if rij + is even (j = 1 or 2), it has a symmetry, which 
is a composition of translation by a half of period along the axis Xj and reflection about 
the mid-plane. The following three combinations of symmetries of this kind have been 
encountered in computations: 

o'l- 7^ only, if mod2ni = mod2n2 = mod2n3 and /in,i 7^ only, if mod2ni ^ 
mod2 712 = mod2 n^; 

C2' Vn,i 7^ only, if mod2ni = mod2n2 = mod2n3 and hn,i 7^ only, if mod2ni = 
mod2 712 7^ mod2 ns; 

(^?>' Vn,i 7^ only, if mod2n2 = mod2n3 and hn,i 7^ only, if mod2n2 7^ mod2n3. (We use 
here the standard notation mod2 n for parity of n: it is for even n, and 1 for odd n.) 

Except for R — 5600 and 5800, the flows in attractors are symmetric about the 
vertical axis (see the table); in agreement with this, magnetic field either has the same 
symmetry, or it is antisymmetric. All the simulated CHM regimes, in which all the three 
constituent fields possess this symmetry, exhibit a simple time dependence: they are 
steady or periodic, and in addition have the symmetries cri, in particular, implying scale 
separation. 

CHM regimes for 5100 < R < 5140 are steady states (see figs. 1 and 2). As documented 
in [25], magnetic field generated by convective fiows tends to concentrate near infinitely 
electrically conducting boundaries of the layer. Magnetic fiux ropes (fig. la) generated by 
steady fiows are associated with stagnation points of the fiow, which have one- dimensional 
unstable and two-dimensional stable manifolds (asymptotic solutions of the magnetic 
induction equation describing magnetic ropes in the vicinity of such stagnation points were 
presented in [15, 29, 8]); the ropes are cut in halves along their axes by the boundaries of 
the layer. In the steady CHM state for R = 5100, stagnation points of the flow are located 
on the midplane and on the upper and lower boundaries of the box of periodicity shown 
on fig. 2. On each of the three horizontal planes they constitute a square mesh of mesh 
size Lh/V^, oriented along the diagonals of the periodicity cells. Only stagnation points 
at the centre and in the vertices of the horizontal faces of the shown box of periodicity 
are associated with the magnetic fiux ropes. The fiow can be described as a system of 
downward and upward swirls (fig. 2) positioned in the chessboard order. Inward fiows 
advect magnetic field from the magnetic ropes near the boundaries towards the middle of 
the layer with formation of magnetic fiux sheets between adjacent swirls (fig. lb). 

A branch of periodic orbits exists for 5160 < R < 5190. (For symmetry reasons, 
this branch of periodic regimes does not emerge from the branch of steady states.) The 
symmetry about the vertical axis of the fiow and temperature is preserved, but magnetic 
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Figure 2: Isosurfaces of kinetic energy density |vp at the level of a half of the maximum 

(a) and the flow on the mid-plane X3 = 0: horizontal components (b) of the velocity of the 
fluid and isolines of the vertical component (c) of the velocity step 0.3 (zero and positive 
values: solid lines, negative values: dashed lines) for the CHM steady state for R = 5100. 
One periodicity cell is shown. 
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Figure 3: Isosurfaces of kinetic energy density |vp at the level of a half of the maximum 
for the time-periodic CHM regime for R = 5180. One periodicity cell at a randomly 
chosen time is shown. 




of the maximum for the time-periodic CHM regime for R = 5180. One periodicity cell is 
shown at the time, for which fig. 3 is constructed. 

field is now antisymmetric (and significantly stronger). Unhke for the CHM steady states 
observed for smaller considered R, there is no parity selection of wave numbers in the 
Fourier series (8)-(10) representing the solutions. However, the structures of the flow 
and (to a lesser extent) magnetic field resembles those for R = 5100 (fig. 3, 4a). The 
main feature of magnetic field is vertical magnetic flux sheets (fig. 4b) , resembling sheets 
developing in the vicinity of an unstable manifold of a stagnation point of the flow, if the 
point has a two-dimensional unstable and one-dimensional stable manifolds (asymptotic 
analysis of solutions of the magnetic induction equation in this context was carried out 
by Childress and Soward 1985). 

At R = 5210 two CHM attractors are found in the system: amagnetic steady rolls and 
an attractor exhibiting heteroclinic chaotic behaviour (fig. 5). In the phases of oscillations 
of growing amplitude the system is in the vicinity of the former periodic orbit. The 
sample trajectory vigorously exponentially departures from it, but subsequently comes to 
the phase of lower kinetic and magnetic energy levels, during which it is in the vicinity 
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Figure 5: Total kinetic (upper curve) and magnetic (lower curve) energies (vertical axis) 
as a function of time (horizontal axis) for the chaotic CHM regime for R — 5210. 

of a mildly unstable steady state, near which integration of the sample trajectory has 
begun. Duration of this phase can be quite short. The phase of evolution near the steady 
state also finishes in an exponential departure, with kinetic and magnetic energy surges 
of amplitude comparable to that of energy surges during departures from the phase of 
evolution near the periodic orbit. Sample trajectories which start in the vicinity of the 
same steady state for higher R, are attracted to steady rolls in shorter times, than the 
length of integration of the trajectory for R — 5210; however, their behaviour at the 
intermediate transitory stage resembles the chaotic behaviour for R = 5210, and thus 
from a conservative point of view it cannot be excluded, that duration of the chaotic 
stage at = 5210 is also finite. 

The window of amagnetic convective rolls is followed by a sequence of periodic regimes 
in the interval 5350 < R < 5550. Periodic CHM regimes from the branch for 5350 < R < 
5550 exhibit an interesting feature: whilst for 5100 < R < 5190 the strongest magnetic 
field is located on boundaries of the convective layer, in the former Rayleigh number 
interval the maximum of magnetic field density is always (5350 < R < 5500) or most of 
the time {R — 5550) is in the interior of the layer (see fig. 6, 7; segments of some plots 
are not visible on fig. 7 for the times, where the distance vanishes). This is in contrast 
with the behaviour found in all previous computations. For < s < 1, denote by Us{t) 
the set of points in the layer, where at time t magnetic energy density |hp is equal to or 
exceeds the value seinax(^), ^maxit) being the maximum of magnetic energy density spatial 
distribution at time t. For R — 5350 and 5400, the set ?7i5/i6(^) is always at a distance 
of at least 0.226 and 0.156 of the layer width away from the boundaries, and Ui/2{t) is 
not less than 0.133 and 0.057 of the width away, respectively. Consequently, these CHM 
regimes may be not very sensitive to the boundary conditions for magnetic field. In all 
simulated CHM regimes the distance is the smallest (in time) near the points of maxima 
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Figure 6: Isosurfaces of magnetic energy density |hp at the level of 1/2 of the maximum 
for the time-periodic CHM regime for it! = 5400. One periodicity cell is shown step 1/8 
of the temporal period T. 
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Figure 6 (the end) 



of the total kinetic and magnetic energies. Concentration of magnetic field in the interior 
of the layer can be related to the fact, that the regions of the most vigorous flow motions 
are inside it (fig. 8, 9). Since isosurfaces shown on these figures only weakly depend on 
time, they arc presented step a half of the temporal period, when the difference in their 
shape is extreme. (Vorticity admits the maximum on the boundaries of the layer; both 
the flow and vorticity are strong in some regions close to the boundaries, but the size 
of these regions is comparatively small, which may explain why they do not influence 
magnetic field generation significantly.) 

The branch existing for 5350 < R < 5550 apparently disappears between R = 5550 
and 5600. For R = 5600 the regime is a periodic orbit of temporal frequency more than 
twice higher (see the table). It has a different symmetry group: it is the only periodic 
orbit that we have found, where the flow is not symmetric about a vertical axis. Between 
R — 5600 and 5700 the periodic regime changes once again. At R — 5800 the CHM regime 
is a chaotic attractor of heteroclinic nature. Like for R = 5210, the sample trajectory 
enjoys excursions between a steady state and a periodic orbit (in this description we 
do not distinguish states related by discrete symmetries and translation in horizontal 
directions), but apart from this the patterns of behaviour of the two chaotic attractors 
are completely different. The periodic orbit emerges from the one for R — 5700 in a 
period-doubling bifurcation; in particular, it has all the parity-related symmetries of the 
latter. This periodic orbit is stable at i? = 5750, but becomes unstable at R = 5800. 
The steady state is amagnetic rolls aligned with the diagonal of the periodicity cell and 
invariant with respect to translations along the axis; therefore, it also possesses all the 
symmetries of the periodic orbit at R — 5700. However, all the symmetries are broken 
during the exponential transition processes when the trajectory moves between the two 
states, except for the composition of translation in the 0:2— direction and reflection about 
the mid-plane a^. Both states are only marginally unstable. The rolls are apparently 
capable of kinematic generation of magnetic field, but while the trajectory approaches 
them, magnetic energy falls off almost to zero (fig. 10). 
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Figure 7: Distances (vertical axis) from sets, where magnetic energy density is equal to 
or exceeds ei„ax(i)/16 (thin short-dashed line), ei„ax(^)/8 (dash-and-dots hne), einax(^)/4 
(dotted line), emax(^)/2 (thin long-dashed line) and ^e^-aa.x{t) (thin solid line) to the closest 
boundaries of the fluid layer for 5350 < R < 5600 as a function of time (horizontal axis). 
Thick solid line: normalised total magnetic energy E"^[t)/ maxtE"^{t), thick dashed line: 
normalised total kinetic energy E^{t)/ maxj E^{t). One temporal period is shown. 
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Figure 8: Isosurfaces of kinetic energy density |vp at the level of a half of the maximum 
for the periodic CHM regime for R — 5400. One periodicity cell is shown step a half of 
the temporal period. 




the periodic CHM regime for R — 5400. One periodicity cell is shown step a half of the 
temporal period. 

4 Conclusion 

We have observed a number of steady and periodic CHM attractors, possessing the symme- 
try about the vertical axis, and several regimes, where this symmetry is broken partially 
(the magnetic field becoming antisymmetric) or completely. The amplitude equations 
derived in [31] are applicable for examination of the weakly non-linear stability of the 
symmetric CHM regimes presented here. It remains to be investigated how disappears 
the branch of periodic orbits, existing for 5350 < R < 5550. If this is a symmetry-breaking 
bifurcation, the system of amphtude equations for free CHM convection, analogous to the 
system of mean-field equations derived in [30] for a branch of regimes emerging in a Hopf 
bifurcation in forced CHM convection, might be applicable for examination of weakly 
nonlinear large-scale stability of the periodic regimes constituting this branch near this 
bifurcation. 

We have found two other types of notable convective magnetic dynamos. One of the 
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Figure 10: Total kinetic (a) and magnetic (b) energies (vertical axes) as functions of time 
(horizontal axis) for the chaotic CHM regime for R — 5800. Magnetic energy is shown in 
logarithmic scale. 
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regimes is periodic orbits for Rayleigh numbers from the interval 5350 < R < 5500: it 
is apparently the first known example of CHM regimes, where the strongest generated 
magnetic field always resides in the interior of the fluid layer. The second is the heterochnic 
chaotic regime for R — 5800, in which the evolution consists of phases, where the regime 
is in the vicinity cither of a periodic orbit, or a steady state - amagnetic rolls. Dynamics 
of this attractor can be studied employing expanded central manifold reduction [18, 19]. 
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